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ABSTRACT 

We have compiled a sample of early-type cluster galaxies from < z < 1.3 and measured the 
evolution of their ellipticity distributions. Our sample contains 487 galaxies in 17 z > 0.3 clusters 
with high-quality space-based imaging and a comparable sample of 210 galaxies in 10 clusters at 
z < 0.05. We select early-type galaxies (elliptical and SO galaxies) that fall within the cluster i?2oo, 
and which lie on the red-sequence in the magnitude range —19.3 > Mb > —21, after correcting 
for luminosity evolution as measured by the fundamental plane. Our ellipticity measurements are 
made in a consistent manner over our whole sample. We perform extensive simulations to quantify 
the systematic and statistical errors, and find that it is crucial to use point-spread function (PSF)- 
corrected model fits; determinations of the ellipticity from Hubble Space Telescope (HST) image data 
that do not account for the PSF "blurring" are systematically and significantly biased to rounder 
ellipticities at redshifts z > 0.3. We find that neither the median ellipticity, nor the shape of the 
ellipticity distribution of cluster early- type galaxies evolves with redshift from z to z > 1 (i.e., 
over the last ^ 8Gyr). The median ellipticity at z > 0.3 is statistically identical with that at z < 0.05, 
being higher by only 0.01 ± 0.02 or 3 ± 6%, while the distribution of ellipticities at z > 0.3 agrees 
with the shape of the z < 0.05 distribution at the 1-2% level (i.e., the probability that they are drawn 
from the same distribution is 98-99%). These results are strongly suggestive of an unchanging overall 
bulge-to-disk ratio distribution for cluster early-type galaxies over the last ~ 8Gyr from z 1 to 
z ^ 0. This result contrasts with that from visual classifications which show that the fraction of 
morphologically-selected disk-dominated early-type galaxies, or SOs, is significantly lower at z > 0.4 
than at z ~ 0. We find that the median disk-dominated early-type, or SO, galaxy has a somewhat 
higher ellipticity at z > 0.3, suggesting that rounder SOs are being assigned as ellipticals. Taking 
the ellipticity measurements and assuming, as in all previous studies, that the intrinsic ellipticity 
distribution of both elliptical and SO galaxies remains constant, then we conclude from the lack of 
evolution in the observed early-type ellipticity distribution that the relative fractions of ellipticals and 
SOs do not evolve from z ~ 1 to z = for a red-sequence selected samples of galaxies in the cores of 
clusters of galaxies. 

Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD, — galaxies: 
evolution — galaxies: fundamental parameters — galaxies: photometry — clusters: 
individual: CL J1226. 9+3332 
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1. INTRODUCTION 

SO galaxies are enigmatic objects whose formation and 
evolution is still not well understood. Originally, these 
galaxies were postulated to exist as a transition cl ass be - 
tween the elliptical and spiral sequence fHu bblelll936f) . 
Both elliptical and SO galaxies lack spiral arms or major 
dust features. However, early work on the ellipticity dis- 
tributions of galaxies showed that SO galaxies were disk- 
dominated systems with ellipticity distributions that dif- 
fered from ellipticals, being more s imilar to spiral galax- 
ies in their intrinsi c shapes {e.g., iRood fc BaumlTl967l : 
ISandage et al.lll970D. The early wor k on the morphology- 
density relation bv lDressle^ (|1980af ) emphasized that SOs 
and ellipticals both occur with higher frequency in higher 
density environme n ts, w hile t he pioneering studies o f 
iButcher fc Oemleil (|1984[ ) and iDressler fc GunnI (119921 ) 
began to provide hints about how the early-type popu- 
lation of elliptical and SO galaxies might evolve out to 
redshifts around 0.5 and earlier. 

GEPI, Observatoire de Paris, Section de Meudon, 5 Place J. 
Janssen, 92195 Meudon Cedex, France; simona.mei@obspm.fr 
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However, it was not until HST allowed comprelien- 
sive hlKh-resolution imaging of distant clusters of galaxies 
that iDressler et all (|1997f ). for example, and others be- 
gan to show directly that the SO fraction was changing 
at high redshift. This early wor k, along with more re- 
cent studies fe.g.. lPostman et a l.l i2005 ,: Desai et al.ll2007f l 
found smaller SO galaxy fractions in clusters of galax- 
ies at higher redshifts, z > 0.3 — 0.4. The implica- 
tion is that the SO galaxy population forms with differ- 
ent time-sc ales and l ater than the ellip ti cal population. 
iFasano et al. (2000' ). ISmith et al.l (l2005l). iPostman et all 
(l2005D,|Poggianti c t all (l2006lfl. iDesai et al.l (|2007t ). anX 
most recently, Wilman et al.l ()2008l ) all find that the ma- 
jority of evolution occurs since z ~ 0.4, i.e., in the past 
4 Gyrs, and above those redshifts there is little or no 
evolution in the early-type galaxy fracti on in clusters of 
galaxies out to z ^ 1 (see ISmith et al.l [2005, for a dif- 
ferent point of view). Since SO and elliptical galaxies 
have different ellipticity distributions and bulge-to-disk 
distributions, two simple (and related) predictions that 
can be drawn from the observed changes in SO fraction 
with redshift is that both the mean bulge-to-disk ratio of 
early-type cluster galaxies and the ellipticity distribution 
should have changed over relatively recent epochs (since 
z ~ 0.4). 

The evidence for evolution in the SO population of clus- 
ters rests primarily on morphological classifications of 
galaxies. The separation of the early-type galaxy popu- 
lation into SO and elliptical galaxies has long been rec- 
ognized as being a challenging task (see, for example, 
lAndrcon 1998 for a discussion of the systematic errors 
in classification and how misclassification mimics evolu- 
tion.) SO galaxies are defined as multicomponent disk 
galaxies, while ellipticals are defined as single compo- 
nent systems. However, a number of quantitative stud- 
ies found that a substantial fraction of ellipticals have 
"disky" isophotes (see iKormendv fc DiorgovskMl989l for 
a sum mary) in contrast wi th this definition. The anal- 
ysis of iRix &: White! (|199Cl[ ) exemplifies the challenge of 
establishing the relative contributions of elliptical and 
SO galaxies. They showed, for a z ^ galaxy, that de- 
tecting a disk component in a spheroidal galaxy is in- 
creasingly difficult as the disk b ecomes more face-on in 
projection. iRix fc White! (|1990( ) found that a disk con- 
taining 20% of the total light of a galaxy is impossible 
to detect over half of the range of cos(i) where i is the 
inclination angle. Specifically, they note "since the cos 
i axis can be interpreted as a probability axis, this im- 
plies that 50% of all disks with Ld/Lb < 0.25 cannot 
be detected by photometric means." These, however, re- 
sults for a single galaxy. For an ensemble of galaxies, 
the bulge-to-disk ratio distribution can be constrained 
by the intrinsic ellipticity distribution, with the average 
observed e llipticity being directly rela ted to the intrinsic 
ellipticity (jBinnev fc Merrifield![l99l ). Therefore, exam- 
ining the distribution of galaxy ellipticities provides a 
direct measure of the evolution in the distribution of the 
b ulge-to-disk ratios o f that g alaxy population. 

!j0rgensen fc Fraip^ (|1994l ) investigated the ellipticity 
distribution of elliptical and SO galaxies in the Coma 
cluster, and suggested that elliptical and SO galaxies 
were not distinct classes but were part of a contin- 
uum of objects of varying bulge-to-disk ratios. This re- 
sult was given additional support by a recent study by 



!Krainovic et^ (|2008D who found that 69% of ellipti- 
cal galaxies have multiple kinematic components, gen- 
erally disk-like components, while 92% of SO galaxies 
have disk-like components. These results give empha- 
sis to the view that elliptical and SO galaxies form a 
continuous distribution of disk frac tion as opposed to 
two d istinct classes. In particular, !j0rgensen fc Fram3 
() 19941 ) constructed a model of the ellipticity distribu- 
tion of elliptical and SO galaxies using a continuum of 
bulg e and disk components vie wed from a variety of an- 
gles. !j0rgensen fc Fraip^ (|1994! ) found a deficit of round 
SO galaxies, in contrast with what was expected from 
their model, suggesting that s ome face-on SO g alaxies had 
been class i fied as ellipticals. iBlakeslee et al.! 120061 and 
iMei et aTl ()2006aD also found a lack of round SO galax- 
ies in three z ~ 1 clusters of galaxies when compared to 
what was expected for a disk population viewed at a va- 
riety of angles. These three studies illustrate a potential 
bias in the visual classification of galaxies, namely that 
nearly face-on SO galaxies galaxies may be incorrectly 
classified as ellipticals. A number of studies suggest that 
quantitative measurements do not suffer the same orien- 
tation bias as visual classifications (IBlakeslee et al.!l2006l : 
Ivan der Wei fc van der Marel! I2008D . but robust, bias- 
free galaxy classification, either visual or quantitative, 
still remains an elusive goal. 

The importance of the work showing an apparent evo- 
lution in the SO fraction in clusters and the knowl e dge o f 
pote ntial classificat i on bia ses led IDressler et al.l ()1997l ) 
and iPostman et al.l (|2005! ) to investigate the ellipticity 
distributions of the SO and elliptical galaxies in their 
samples. In general, the authors found that the ellip- 
ticity distributions of SO and elliptical galaxies show no 
evolution over the bro ad redshift ranges in their samples 
(jPostman et al.! [20051 compared clusters from z ^ 0.25 
to z ~ 1.3). Also, the ellipticity distributions of elliptical 
and SO galaxies differ from each other, providing evi- 
dence for the existence of two distinct classes of galaxies. 
However, these previous studies do not use a consistent 
measure of the ellipticity as compared wi th the z ^ ef- 
forts s uch as[j0rgensen fc Franx (199 J) or lAndreon et al.! 
(1996!) , making a comparison between these z > 0.2 
samples observed with HST and z ~ samples ob- 
served from the ground difficult. Traditionally, elliptici- 
ties for galaxies at z ~ were measured by visual esti- 
mates (Dressier 19 80al) or by fitting models to the ellip- 
tical isophotes (jJorgensen fc Franxlll994HAndreon et al.! 



1996t). At higher reds hifts, z > 0.2, t he elh pticities in 



Small et all (|1997f) or IPostman et al.! (|2005D are deter 



mined by the second-order flux-weighted moments of a 
detection isophote. 

Both the data and the techniques have matured so that 
we can now evaluate the ellipticity distribution, as a func- 
tion of redshift, quantitatively, and even more impor- 
tantly, in a consistent way with minimal systematic er- 
ror. We show that some approaches used previously that 
could not correct for the point-spread function (PSF) are 
probably subject to significant systematic error. Essen- 
tially, the "blurring" effect of the PSF will lead to galax- 
ies being measured as rounder than they actually are. We 
will use a single consistent and robust approach for mea- 
suring the ellipticities at all redshifts. The ellipticities 
will be measured by fitting models convolved with the 
PSF to galaxy surface brightness profiles. This will elimi- 
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nate some of the previously-reported uncertainties found 
when comparing ground-based im aging data taken un- 
der different seeing conditions fsee lAndreon et al]|1996l 
for some discussion). In addition, this approach essen- 
tially eliminates the systematic error associated with the 
PSF "blurring" . Our ellipticity measurements will pro- 
vide an assessment of the evolution in the distribution of 
the overall bulge-to-disk ratio of early-type cluster galax- 
ies from the present day to redshifts z ^ 1. 

We will discuss how we compiled our samples of early- 
type cluster galaxies with morphological classifications in 
Sj2l and then how we measured their total magnitudes, 
colors, and ellipticities in fj3l One of the advantages of 
our approach of using an automated measurement tech- 
nique is that we can simulate the measurement process. 
We discuss this in !j3.3l and in Appendix Al. From our 
measurements, we find no evolution in the distribution of 
ellipticities of cluster early-type galaxies, which we show 
in SjH In fjsl we discuss the implications of this result 
and contrast with previous measurements of the evolu- 
tion in the overall distribution of the bulge-to-disk ratio 
of cluster early-type galaxy population. We follow this 
with a summary and some discussion of the broader im- 
plications in Sjni 

Throughout this paper, we assume fl„i = 0.27, ^Ia — 
0.73 and Ho = 71 km s~^ Mpc~^. All B magnitudes we 
list use the Vega zeropoint {Bvega — + 0.11.) Other 
magnitudes where given are AB mags. 

2. SAMPLE SELECTION 

To carry out our study of the ellipticities of early-type 
galaxies from z 1 to the current epoch, we have assem- 
bled a sample of morphologically-selected, cluster early- 
type galaxies ranging in redshifts from z ~ to 1.27. 
The clusters used are tabulated below, as is the source 
of the morphologies for the galaxies samples ( g2.2p . The 
early-type galaxies are chosen to lie within the projected 
R200, which we use to define cluster membership. The 
measurement for each galaxy of its color, magnitude, and 
ellipticity, is described in ^ 

2.1. Cluster Sample Selection 

Our z < 0.05 sample is selected from Abell clus- 
ters in the Sloan Digital Sk y Survey Fifth Data Release 
(|Adelman-McCarthv et alH 2007. SDSS DR5). We sum- 
marize t his sample in Tab le [TJ where we list the Abell 
number ( Abell et al.lll989| ) and the redshift of the clus- 
ter from Struble &: RoodI (fl999l) . The last column lists 
the number of redshift-selected members with early-type 
classification, which lie within the projected i?200 and are 
on the g — r red sequence in the SDSS DR5 imaging. 

To determin e i?20C )i we will use the formula give in 
ICarlberg et~al\ (|1997f l. or 



-R200 — 



V3 



CTi 



10 H{z) 



where ui is the one-dimensional velocity dispersion and 
H{z) is th e Hubble c o nstant at the redshift of ob- 
servation. iDesai et al.l (|2007t ) uses the relation from 
iFinn et all (120051 ) . which is functionally the same, so ef- 
fectively our selection radii are similar to iDesai et al.l 
(|2b07.) . We selected 2i?20o/7'', instead of i?200, as a 
galaxy at -R200 from the cluster center will be, on av- 



erage, projected to appea r at the distance 2R2qo/'k (see 
iLimber fc Mathewill960l . for example). 

For the high redshift clusters, we required Hubble Space 
Telescope (HST) imaging with the Advanced Camera for 
Surveys (ACS) or the Wide Field Planetary Camera 2 
(WFPC2) in multiple filters. We also required that the 
morphological classificati ons were done i n a way con- 
sistent with the original iDressleil (|1980bD and iDressleil 
(19803) work. For the three clusters where 2R2qo/it 
was larger than the field of view over which we had 
imaging data (gene rally the clusters in the sample of 
iDressler et al.lll997l ). we simply used all available galax- 
ies. We mark those clusters in Table [2] with the super- 
script b. Below we will detail the sources for the mor- 
phological classifications and how we measured the total 
magnitudes, colors, and ellipticities. 

We tabulate the 0.3 < z < 1.3 clusters for which we 
have early- type galaxy samples in Tabled In that table, 
we list the dispersions, inferred radii, and the final sample 
sizes. For CL J1226-1-33 we have compiled a new catalog 
of members which we will discuss in Appendix [Bl For 
RX J0849+4452, we compute a ne w dispersion, ^ye use 
both redshifts fr om previous work (jStanford et al.ll2001l : 
iMei et al.ll2006bf) and unpublished ones we have recently 
collected (which we plan to publish in a future paper.) 
There are a total of 18 galaxies in RX J0849-I-4452 with 
redshifts within i?200- The biweight center of the distri- 
bution is z = 1.2600 ± 0.0017 and the dispersion is 798 
± 208 km s^^. The errors for both the redshift and the 
dispersion are estimated by bootstrapping the redshift 
distribution. 

For the clusters RX J0152-13, MS 1054-03, MS 2053- 
04, and CL 1358- 1- 62, w e used the redshift catalogs 
from iHolden et all (|2007t ). The faint magnitude limit 
we adopted corresponds to the completeness limits for 
our high redshift samples (Holdcn et al. 2007). For the 
remaining clusters, we included all galaxies that met our 
magnitude, color, morphology and radial requirements, 
rejecting only those galaxies with redshifts outside of the 
clusters. This will mean that some cluster samples will be 
contaminated, i.e., some red-sequence early-type galax- 
ies will not be members. The level of contamination is 
discussed below. Our z > 0.3 sample consists of 487 
early-type galaxies while our z < 0.05 sample has 210 
galaxies. 

Many of the cluster red-sequence selections we use are 
tabulated in other studies. For those clusters studied by 
the ACS Instrument Definit ion Team, these red-sequence 
relations are summarized in lMei et al.l ()2008D . For the re- 
maining clusters, we derive the early- type red-sequence 
from the data. Whenever possible, we use existing red- 
shift catalogs to determine the red-sequence of spectro- 
scopic determined members. We then accept all galaxies 
that lie within the 2fj of that sequence, rejecting those 
galaxies that are known not to be members. 

We know that there will be some contamination by 
field galaxies in our catalogs. To measure the level of con- 
tamination, we used our redshift catalogs for MS 1054- 
03 ai id RX J0152-13 along with the catalogs of EDisCS 
from iHallidav et al.l (|2004l i. RX J0152-13 has a known 
group in the foregroun d with colors very similar to 
that of the cluster (see IHolden et al] [20f)5l for a dis- 
cussion). Nonetheless, when we consider all early- type 
galaxies that lie on the red-sequence for MS 1054-03 and 
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TABLE 1 

z < 0.05 Cluster Summary 



Cluster z cr^ -R200'' # Early-type'^ 

(kms-l) (Mpc) 



ACQ 119 


0.0440 


744 


1.13 


41 


ACQ 168 


0.0452 


524 


0.80 


54 


ACQ 194 


0.0178 


435 


0.67 


9 


ACQ 957 


0.0440 


691 


1.05 


32 


ACQ 1139 


0.0383 


436<* 


0.67 


23 


ACQ 1142 


0.0353 


417 


0.64 


17 


ACQ 1656 


0.023 


1008 


1.55 


203 


ACQ 1983 


0.0441 


433 


0.66 


18 


ACQ 2040 


0.0456 


602 


0.92 


50 


ACQ 2063 


0.0337 


521 


0.80 


50 


ACQ 2151 


0.0371 


786 


1.20 


58 



The cluster redshift and a from IStruble fc RoodI (119991) 
unless otherwise noted. ^ R200 is derived from The 
number of galaxies with E, E/SO, SO/E, SO, or S/SO 
classifications.'* From iPoegianti et al.l II2006I ) 



RX J0152-13, regardless of redshift, we find a 3 ± 1% 
contamination rate of nonmembers. When we examine 
the less extremely rich sample of the EDisCS, we find a 
contamination fraction of 10 ± 3%, which is likely more 
representative of the typical clusters in our sample. The 
small size of this contamination means we need not com- 
pute a statistical background correction for the cluster 
red-sequence. 

2.2. Morphologies 

Galaxy morphologies were obtained from the literature 
for all th e galaxies in our sample. We use the morpholo- 
gies from lDressleij (|1980al ) for the z < 0.05 cluster mem- 
bers. We removed all late-type galaxies from that low 
redshift sample, i.e., all those classified as spiral, irreg- 
ular or unknown classifications. We considered galaxies 
classified as SO/a as SO, and galaxies classified as Sa/0 as 
spirals. We list in Table[l]the number of early-type galax- 
ies - E through SO/a - in the last colu mn. The numbers 
are smaller than the number available in lDressleil (|1980af) 
because not all of the galaxies in each cluster are imaged 
in the SDSS DR5. 

For o ur higher red shift clusters, we used t he classifica- 



Dressler et al. (^19973 as tab ulated inlSmail et al.l 

' ool 



Postman et all (|2005l ). and iDesai et all (|2007f ). 



tions of 
([1993), 

Each cluster was imaged with either ACS or WFPC2 
with morphological classifications d one in a. manne r 
consistent with the previous work of iDressleH (|1980aD . 
The clusters MS 1054-03, RX J0152-13, CL J1226-F33, 
CL 1604-1-4304, CL 1604-1-4321, RX J0910-I-5422, 
RX J1252-292 7, and RX J0849-f4452 all come from 
iPostman et"an ()2005[) . The MORPHS survey clus- 
ter s CL 0016-H6 , CL 054-27, CL 0024-1-16 are part 
of iDressler et al.l ()1997[ ) sample. Five clusters from 
the EDisCS survey - ECL J1040-11, ECL J1054-11, 
ECL J1054-12, ECL J 1216-12 a nd EC L J1232-12- have 
morphologies from De sai et al.l (|2007( ). Two additional, 
CL 1358-F62 a n d MS 2053 -04, have morphologies in 
iFabricant et~al] (|2000l ) and iTran et al.l (2001), respec- 
tively. 

3. DATA AND MEASUREMENTS 

Our sample of early-type galaxies is based on a range 
in luminosity, corrected for the observed passive evolu- 



tion to z ~ 1. The derivation of the total magnitudes 
and colors used for our sample selection is discussed be- 
low. Taking that luminosity-selected sample, we then 
discuss the derivation of their ellipticities. These con- 
stitute the key observable for this work, and the impli- 
cations of those ellipticity measurements are discussed 
in the remainder of the paper. Our final selection will 
consist of morphologically-selected early-type galaxies on 
the red-sequence within a well-defined magnitude range 
(—19.3 > Mb > —21 at z = 0) lying within the cluster 
core that is defined by i?200 ■ These magnitude limits cor- 
respond to Mq + 1 > Mb > - 0.75 and to galaxies 
with stellar masses roughly between 10^°'^ Mq < M < 
10^^-^ Mq, assu ming a "diet" Sa lpeter initial mass func- 
tion (IMF) (see iBeU et aLll2003l for a discussion of the 
IMF and the procedure we use to estimate the stellar 
masses of galaxies). 

3.1. Total Magnitudes and Colors 

Our samples were selected based on the rest-frame B 
magnitude. To estimate these magnitudes, we used the 
total magnitude in the passband closest to the rest-frame 
B. We also needed a color to correct the apparent mag- 
nitude in the observed passband to a rest-frame B mag- 
nitude; how this was done is outlined below. The Sersic 
model fits that we used to measure the ellipticities also 
were used to determine the total magnitudes and the 
color apertures. The total magnitude is the normaliza- 
tion of the Sersic model fit. For the color aperture, we 
used the circularized half-light radius, rhir — ahir^/Q, 
where q is the ratio of the minor to major axis, or 1 — e, 
and ahir is the half-light radius along the major axis of 
the best-fitting elliptical model as determine d by GAL- 
FIT. These are the same apertures used in iMei et aTl 
(|2008f) (see that paper for more detail). 

We adjust the magnitudes of the galaxies by 1.208z as 
measu red for early-type galaxies using the fundamental 
plane (jvan Dokkum fc van der Mare]ll200l . This com- 
pensates for the mean passive evolution of the old stellar 
population. Our magnitude range covers M* — 0.7 to 
M* + 1 using the M* froi n iNorberg et a l. (2002) after 
converting the bj used by INorberg et al.l (.2002) to the 
B of the Johnson-Morgan system ( Buseilll978l ) that we 
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TABLE 2 

Summary of 0.3 < 2 < 1.3 Cluster Data 



Cluster 


z 


R.A. 


Dec. 


Obs. Filters 


o- -R200 # Early-type 












(km s^i) (Mpc) 



CL 1358-f 62 


0.328^= 


13 59 50.6 


+62 59 05 


Veoe 


-^814 


1027+45° 


1.4 


67 


CL 00244-16 


0.395 


00 26 35.7 


+17 09 46 


^•625 »775 


2850 




0.83'' 


26 


CL 0016+16 


0.541 


00 18 33.5 


+16 26 14 


«775 


2850 


1234^1^8 d 


1.46*^ 


42 


ECL J1232-12 


0.541° 


12 32 30.3 


-12 50 36 


V I, 


-^814 


1080^99^° 


1.28 


21 


CL 0054-27 


0.56 


00 56 56.9 


-27 40 30 


Veo6 


^814 


1180*^ 


1.38*^ 


19 


MS 2053-04 


0.587 


20 56 21.3 


-04 37 51 


Veo6 


^814 




1.00 


31 


ECL J1054-11 


0.697'' 


10 54 24.4 


-11 46 19 


V I, 


^814 




0.64 


5 


ECL J1040-11 


0.704== 


10 40 40.3 


-11 56 04 


V I, 


-^814 


418^45 ° 


0.45 


2 


ECL J1054-12 


0.750^= 


10 54 43.5 


-12 45 52 


V I, 


-^^814 


bOAtll' ° 


0.53 


13 


ECL J1216-12 


0.794° 


12 16 45.3 


-12 01 18 


V I, 


-^814 


1080t^^*' ° 


1.04 


31 


MS 1054-03 
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Fig. 1. — EUipticity measured from the simulations, which used real galaxy images, as a function of the ellipticity in the original input 
image. The error bars show the scatter around the recovered ellipticity values for each galaxy. Each galaxy image is rescaled in size 
appropriately for the redshift of each cluster. Each galaxy is realized at a variety of signal-to-noise values covering the range over which 
the galaxies in the high redshift sample are observed, and were then placed in the cluster imaging data. We plot a straight line with a 
slope of one, the expected relation if we recover the input ellipticities. The scatter is typically ere ~ 0.01 — 0.03. The scatter increase at 
lower signal-to-noise, and matches the statistical errors from the fitting process. The median offset is —0.01 with a range of 0.00 to —0.03, 
showing that our method accurately recovers the ellipticity of these low-redshift galaxies when they are observed at high redshift. No 
systematic trends are seen as a function of ellipticity or redshift. 
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use in this paper. We trim our sample at Mb < —19.3, 
as our high redshift samples become incomplete fainter 
than that magnitude. The brighter magnitude limit is set 
to be Mb < —21 to exclude the most luminous galax- 
ies (whose formation and evolution may differ, as is ex- 
plained later in M.ip . Thus our adopted magnitude range 
is -19.3 > Mb - 1.208z > -21. 

Using our simulations of real galaxies that we discuss 
in more detail in Appendix IA.21 we estimate the typical 
error and offset for these total magnitude measurements. 
We find that total magnitudes as measured with ACS 
have an error of cr = 0.10 mag while magnitudes mea- 
sured with WFPC2 have a — 0.15 mag. There is a small 
offset at most redshifts, such that we measure a magni- 
tude brighter than the actual magnitude of the simulated 
galaxy. For the clusters at 0.3 < z < 0.6, this is only 0.02 
mag, while it increases to 0.06 mag at 0.6 < z < 1 and 
0.10 mag at 2; > 1. This offset is the same for galaxies 
regardless of size, morphological type or the value of n 
from the Sersic model fit. We do not apply this offset to 
the measured magnitudes. 

To ensure a reliable measure of the c olor, we applied 
the CLEAN algorithm (jHogbomI Il974f ) to the original 
images in all passbands. Using CLEAN mitigates the ef- 
fects of the different- sized PSF in different passbands, see 
iSirianni et al.l (|2005f ) for examples involving ACS. The 
final "CLEANed" images were used for measuring the 
galaxy colors within the r^ir given above. 

3.2. Redshifted Magnitudes 

We transformed the observed magnitudes into 
redshifted magnitude s usin g the same pro c ess a s 
van Dokkum fc Franxl (1 19961) . iBlakeslee et all (|2006D . 



Holden et all (|2006[ ) and iHolden et al.l (|2007D . We cal 



culated the magnitudes of templates in the rest-frame 
filters. We then redshifted the templates, and computed 
the magnitudes in the observed filters. For the tem- 
plates, we used exponentially decay i ng st ar-formation 
rate models from lBruzual fc Charloti (I2003I. BC OS): the 
same models were used in IHolden et al.l ()2OO70 . These 
models had exponential time-scales of 0.1- 5 Gyr, cov- 
ering a range of ages from 0.5 Gyr to 12 Gyr and three 
metal abundances, 2.5 solar, 1.0 solar and 0.4 solar. For 
the re s t-fram e filters, we used the B and V curves from 
iBuseii ()1978 ^). specifically the B3 curve for the Bz as tab- 
ulated by BC03. We use the same templates and proce- 
dure for all of the clusters in our sample. 

3.3. EUipticities 

We measured the ellipticitics for our galaxies using the 
resuhs from GALFIT (Peng et al. 2002). GALFIT fits 
an elliptical model to the surface brightness profile. Ef- 
fectively, the ellipticity measurements we use are ellip- 
ticitics at the half-light radius. The model is convolved 
with a PSF before it is compared with the data. The ad- 
vantages of the approach we have used is that the GAL- 
FIT fit procedure minimizes the effect of the smoothing 
from the PSF. PSF "blurring" will make galaxies ap- 
pear rounder than they actually are, unless the galaxies 
being fitted have sizes much greater than the PSF. In 
Appendix 1X1 we discuss in more detail the robustness of 
the ellipticity measurements. 

This procedure is performed for galaxies at all red- 
shifts, but requires an estimate of the PSF. For the clus- 



ters of galaxies at z < 0.05, our imaging data came 
SDSS DR5. We use the software tools provided by the 
SDSS to extract the PSF appropriate for each galaxy. 
For the higher redshift galaxies, the pipeline processing 
system we use for ACS (Apsis - ACS pipeline science 
investigation software ) prov ided us with suitable PSF 
(jBlakeslee et al.l l2003l I2006D . We used empirical PSF 
models constructed from multiple ACS observations of 
47 Tuc. For the WFPC2 data, the PSF stiU has a strong 
positional dependence and is significantly under sampled. 
For these data, we estimated the PSF using the TinyTim 
software package ()Kristj|1995[) for each galaxy we fit. We 
also used the option in GALFIT to convolve the model 
of the galaxy with an over-sampled PSF. Before compar- 
ing the PSF convolved model with the data, GALFIT re- 
binned the model to the WFPC pixel scale and smoothed 
the model with the charge diffusion kernel. 

3.3.1. Simulations of Ellipticity Measurements 

We constructed a set of simulations using real galaxy 
images and placed those images in the data of the high 
redshift clusters in our sample. This is discussed in detail 
in Appendix IA.2I 

To summarize, we made multiple measurements of each 
simulated individual galaxy to assess both the impact of 
noise and any systematic effects, and each galaxy was 
simulated over a range in magnitudes. We found that 
the scatter in the ellipticity measurements of the im- 
ages of real galaxies was (Te ~ 0.01 — 0.03 at magnitudes 
typical of those in our samples, and only increased (to 
CTe ~ 0.05 — 0.06) at or below the magnitude hmit of our 
samples. These estimates of the uncertainty are in good 
agreement with the errors estimated by the model fitting 
process used to measure the ellipticity. We found no 
large systematic trends in the ellipticity measurements. 
This can be see in Figure [1] where we show the recovered 
ellipticity as a function of the input ellipticity for three 
redshift bins. There was only a very small overall shift 
of 6e — 0.01, or 3% for the typical galaxy, such that the 
typical z > 0.3 galaxy is measured to be slightly rounder 
than it would appear in the z < 0.05 sample. We found 
the shift was the same, regardless of the redshift of the 
galaxy. No systematic changes were seen as a function of 
ellipticity. For the rest of the paper, we ignore the small 
systematic shift of 6e — 0.01 (noted above and in Figure 
1), and simply quote the observed values. 

We also tested the robustness of the ellipticity mea- 
surements by using incorrect PSFs. We found that this 
error could result in a systematic offset of — ±0.01 
to Se = ±0.03 depending on how bad the mismatch be- 
tween the PSFs were. The small size of our statistical 
errors means that we are sensitive to systematic errors 
of this order, and we will discuss the implications of this 
in later sections. 

4. ELLIPTICITY DISTRIBUTIONS 
4.1. Ellipticity Distributions for the z < 0.05 Sample 

In this section, we define our selection of z < 0.05 
galaxies, and discuss how we will characterize the ellip- 
ticity distributions in order to compare them with the 
z > 0.3 cluster galaxy sample. 

In Figure [H we plot the ellipticity as a function of 
Mb, and also by stellar mass, for all galaxies in the z < 
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Fig. 2. — EUipticity versus absolute B magnitude or stellar 
mass for the z < 0.05 clusters, as described in the text in [[3] 
The dashed lines show magnitude limits we will use for this pa- 
per, -19.3 > Mb > -21. Below Mb = -19.3 our higher red- 
shift cluster samples become incomplete, while above brighter than 
Mb ~ —21, the ellipticity distribution of the early-type galaxy 
population changes, becoming rounder, possibly the result of the 
somewhat different evolutionary history of the most massive galax- 
ies. The galaxies within our selection limits are filled circles, while 
all of the remaining data are shown as open circles. The median el- 
lipticity of the galaxies in our sample is shown by a solid line. The 
mass-selected sample and the luminosity-selected samples yield a 
very similar selection with similar median ellipticity, suggesting 
that mass-dependent effects are likely to be small. For the rest of 
the paper, we will use the B magnitude limits in the left panel to 
select galaxies at higher redshifts. 

0.05 sample that are within 2i?20o/i'j a-nd classified as 
an early-type galaxy. We show the median ellipticity of 
the whole population with a solid line. The stellar masse s 
are derived using the prescription from lBell et al.l (|2003f ). 
which uses the rest-frame colors to estimate the mass-to- 
light ratio of the stellar population, assuming a "diet" 
Salpeter IMF. The "diet" Salpeter IMF is a Salpeter IMF 
with truncation at very low masses, resulting in a mass- 
to-light ratio of the stellar population 0.15 dex smaller 
than a Salpeter IMF. These stellar mass estimates agree 
well with the mass estim ates from the fundamental plane 
(see iHolden et aDl2007t ) . 

At the highest masses or brightest magnitudes, the 
pop ulation becomes ro u nder, a result seen in oth er work 
(see iFranx et al.l 119911 : IVincent fc RvdenI l2005l . for ex- 
ample), likely a result of a different morphological mix 
among the most luminous galaxies which may have a 
somewhat different formation and evolution history. We 
will exclude these most luminous galaxies, Mb < —21, 
from our sample. We illustrate our sample range —19.3 > 
Mb > — 21 with dashed lines in Figure [2] 

At low redshifts, we know that elliptical and SO galax- 
ies have different ellipticity distributions. In Figure[31 we 
show that we reproduce those different distributions with 
our z < 0.05 cluster galaxy sample. The elliptical pop- 
ulation, shown in red in Figure [Sj is much rounder than 
the SO population, shown in blue. The whole of the pop- 
ulation is shown in green. Because SO galaxies dominate 
the population, the green line appears closer to the SO 
distribution than the elliptical distribution. The median 
ellipticities are different, with the median ellipticity of SO 



galaxies Cmed = 0.38 ± 0.02 while the median ellipticity 
of elhptical galaxies is emed = 0.18 ± 0.01. 

In Figureini there is a hint of a deficit of round S O galax- 
ies, a result found bv lJ0rgensen fc Franxl ()1994) . We fit 
to the distribution of SO ellip ticities with th e disk galaxy 
model used in J0rg ensen fc Franxl (jl994l ). This disk 
galaxy distribution is that of an oblate spheroid with the 
minor to major axis ratios drawn from a Gaussian distri- 
bution with mean e and a dispersion of . We find that 
the values for the G aussian that best describe s the data 
are similar to that of |j0rgensen fc Franxl ()1994l ). Our val- 
ues of e = 0.73 ± 0.10 and = 0.10 ±0.02 are consistent 
with £ = 0.65 and = 0.10 from iJorgensen fc Franxl 
(Il99l - a differ ence ~ la. At small ellipticity values, 
the ellipticity distribution of visually classified SO galax- 
ies fall consistently below the expectation of the model, 
which indicates a possible lack of round SO galaxies. 
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Fig. 3. — EUipticity distributions for magnitude-selected samples 
by galaxy type for the z < 0.05 clusters. We show the elliptical and 
SO population separately (red - ellipticals, blue - SO galaxies, green 
- all early-type galaxies). The two sets of galaxies show different 
ellipticity distributions in our data, as expected from the results 
of previous work. The elliptical galaxy population shows a peaked 
distribution, with a median ellipticity of e^e^j = O.lSitO.Ol. In con- 
trast, the SO galaxies have a broader ellipticity distribution, with 
^med = 0-38 ± 0.02, consistent with a more disk dominated popu- 
lation. We show, with a black line, the best-fitting disk population 
drawn from a Gaussian distribution which has a mean thickness of 
b = 0.27 it 0.10 giving a mean ellipticity of e = 0.73 ± 0.10, with a 
standard deviation of CTc = 0.10 it 0.02. If the SO fraction decreases 
with redshift, the overall ellipticity distribution of the cluster pop- 
ulation should also evolve as there will be fewer galaxies with large 
ellipticities.. 



Coma (A1656) is a dominant contributor to the early- 
type sample at z < 0.05, with about 37% of the total. 
We investigated the impact of the removal of Coma from 
the sample. As can be seen in Table [3] doing so makes 
surprisingly little difference to the ellipticity measures 
of the early-types or of the E and SOs separately. The 
ellipticity is consistent within the uncertainties. The SO 
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fraction of the z < 0.05 sample is 68 ± 3% for the whole 
sample and 72 ± 4% withou t Coma . 

Bot h |j0rgensen &: Fran;^ ()1994D and lAndreon et all 
(|1996f ) tabulate their ellip ticities. The media i i ellip tici- 
ties in the Coma sample of lJeirgensen fc Franxl (|1994l ) are 
Cmed = 0.34 ± 0.02 for the SO galaxies and €med = 0.16 ± 
0.02 for ellipticals. Both medians are slightly rounder 
than our measured low redshift values, though the dif- 
ferences are small. The lower values are possibly due to 
the use of the PSF in our analysis. iJorgensen fc Fram3 
(fl994i) did not remove the smoothing of the PSF, and so 
w ill be offset to slightly rounder values as a result. 

lAndreon et al.l (|1996l ) derive ellipticies for the SO galax- 
ies and ellipticals. These values have been used as a low- 
red shift c omparison s et by other higher redshift stud- 
ies ([Small e t al. 19971 iFasano et al.|[2000f) . it is useful to 
understand the difference between these measurements 
and ours. These ellipticities are measured at a fixed 
jifi isophote. Thus m easurement of the ellipticity from 
lAndreon et all (1996) is not directly comparable to ours 
and is not the optimal approach for our study. 

4.2. z > 0.3 Ellipticity Distributions 

We plot in Figure [3] the distribution of ellipticities 
as a function of absolute B magnitude for four redshift 
bins. As before, we select only those galaxies in a fixed 
magnitude range. We adjust the magnitudes, however, 
by I.2O82: mag to reflect the amount of passive galaxy 
evolution that is measured from the fundamental plane 
(jvan Dokkum &: van der Marelll20^7f l. We combine each 
cluster's sample into three redshift bins, 0.3 < z < 0.6, 
0.6 < z < 1.0 and z > 1. We show, with a solid line, the 
median ellipticity for the selected sample in each bin. 

We find no evolution in the median ellipticity with red- 
shift. In Figure [5l we plot the median ellipticity for each 
cluster. In this plot we show the median ellipticity for 
each cluster as open circles. For our z < 0.05 sam- 
ple, we show just the median value (the blue square), 
Cmed = 0.29 ± 0.02. The median values for the individ- 
ual z > 0.3 clusters show larger scatter, but in general 
are quite consistent with the low redshift value, with 
the whole z > 0.3 sample having emed = 0.30 ± 0.01. 
The green, orange, and red points are the median ellip- 
ticities (shown also in Figure 2] as the horizontal lines) 
at 0.3 < z < 0.6 and at 0.6 < z < 1.0 and z > 1, 

emed = 0.31±0.02, eraed = 0.30±0.02, emed = 0.29±0.03, 

respectively. The errors we quote, here and later in the 
paper, are the errors on the median, or •\/(7r/2) the er- 
ror on the mean. We confirmed these error estimates 
with bootstrap resampling. We compared the scatter in 
the high redshift sample by computing the around 
the low redshift median value. We find a, xt — 1-27 for 
V = 17 degrees of freedom, confirming both the good vi- 
sual agreement between the high redshift data and that 
the scatter is not higher than expected from random er- 
rors. As we discuss in Appendix IA.2i we expect the 
largest systematic errors on the median ellipticity val- 
ues to be ~ 0.02, which are not large enough to shift the 
high redshift data to significantly smaller ellipticity val- 
ues. The lack of any change in ellipticity with redshift is 
striking. Formally, the median ellipticity of our sample 
of 4S7 early-type galaxies at z > 0.3 is statistically iden- 
tical to that of the 210 early-type galaxies at z < 0.05, 
being higher by only 0.01 ± 0.02 or 3 ± 6%. 



To add a more detailed assessment of the changes and 
to provide a more quantitative basis for the results seen 
above, we plot the ellipticity distributions of the sam- 
ples in Figure [6l both differentially and in cumulative 
form. In each figure, the blue line shows the z < 0.05 
comparison sample, with the normalization rescaled to 
match the higher redshift samples. At no redshift do 
we find a statistically significant change in the elliptic- 
ity distribution. We use a number of tests to quantify 
this, including a Kolmogorov-Smirnov, Wilcoxon-Mann- 
Whitney rank sum test and a Kuiper test. Taking the 
results from these tests we find that the distribution of el- 
lipticities at z > 0.3 agrees with the shape of the z < 0.05 
distribution at the 1-2% level ( i.e., the probability that 
they are drawn from the same distribution is 98-99%). 
This is a remarkable demonstration of the consistency of 
the ellipticity distributions over a time span of more than 
half the age of the universe. 

There is a hint, in our z > 1 sample, of a lack of high 
ellipticity galaxies, though it is not significant in any of 
our tests. This is, in large part, because of the small sam- 
ple size. There are only 48 galaxies in the z > 1 cluster 
sample. In M.l) we showed that the median ellipticity 
for the z < 0.05 sample was essentially unchanged when 
Coma (A1656) was removed (see Table [3]). As expected 
from that test, there is essentially no change in these re- 
sults on the ellipticity trends whether Coma is included 
or removed from the z < 0.05 sample. 

5. DISCUSSION 

We find that neither the median ellipticity nor the 
shape of the high redshift ellipticity distribution of early- 
type cluster galaxies evolves with redshift, implying no 
change in the overall distribution of the bulge-to-disk ra- 
tio of early- type galaxies with redshift. As we show in 
FigurelSl the SO population at low redshift (z < 0.05) has 
a different ellipticity distribution than the elliptical pop- 
ulation, with SO galaxies having a higher median elliptic- 
ity. If the SO fraction of the early-type galaxy population 
decreases, then the median ellipticity of the early-type 
population should decrease. The lack of evolution we 
observe in the median ellipticity and in the shape of the 
ellipticity distribution implies little or no evolution in the 
SO fraction. This differs fro m the expecta. t ion fr om pre- 
vious work, such as that of iDressler et al.l (|1997f ). which 
finds a decrease in the SO fraction with redshift. 

5.1. Morphological Evolution in the z > 0.3 Cluster 
Sample 

One possible reason why we may not find any evolu- 
tion in the ellipticity distribution of the cluster early-type 
population could be because of the nature of our sam- 
ple. The fraction of galaxies morphologically selected as 
SO galaxies may not change in our sample as has been 
found by other authors. However, when we look at the 
SO fraction of our sample with redshift, using the vi- 
sual classifications for our sample from the literature, we 
find a similar trend to what has been reported in other 
papers, i.e., a lower fraction of SO galaxies at redshifts 
z > 0.4. In Figure [71 we plot the fraction of morphologi- 
cally identified SO galaxies with redshift (from the studies 
discussed in ^i^. O ur sample has proper t ies consistent 
with previous work (IDressler et al.lll997t iFasano et al.l 



l2000tlPostman et al.ll2005tlDesai et al.ll20071 )" as would be 
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TABLE 3 

Median Ellipticities of 2 < 0.05 Sample 



Sample 


All Early-types 


SO's 


Ellipticals 


All z < 0.05 clusters 


0.29 ± 0.02 


0.38 ± 0.02 


0.18 ± 0.01 


Without Coma 


0.29 ± 0.02 


0.34 ± 0.03 


0.21 ± 0.02 
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Fig. 4. — EUipticity versus absolute B magnitude for all the clusters in our sample. The magnitude range, shown by the dashed lines, is the 
same throughout, as we derived in Figure[2] Each galaxy is selected to be within 2/2200/"" ^^'^ to have —19. 3 > Mg + 1.208^ > —21, after 
remo ving the effects of passive evolution. We assume early-type galaxies become brighter by 1.2082 mag llvan Dokkum &: van der Marell 
1 20071 V In each panel, the solid dots are those in the magnitude range of our selection, while the open circles show the remaining galaxies 
in our sample. Each panel covers a different range in redshifts, the leftmost, z < 0.05, followed by 0.3 < z < 0.6, 0.6 < z < 1.0, and 
2 > 1.0. The horizontal line is the estimate of the median ellipticity for the sample within the magnitude limits. It is striking to note that 
the median ellipticity distribution does not change with redshift. Only at the highest redshifts do we see a hint of fewer highly elliptical 
galaxies. 



expected since our sample largely overlaps with previous 
studies and uses visual classifications from those studies. 
We find 42 ± 2% of the early-type galaxy population are 
SO galaxies at z > 0.4. For comp arison, the frac t ion o f 
SO galaxies within 2i?2oo/''' for the lPostman et all ()2005f ) 
"z ^ 1 composite" sample is 35 ± 3%. This shows that 
our red-sequence selection is consistent, being less than 
2(7 different (the Postman et al. 2005 value is derived by 
summing the individual listings in that paper's Table 4 
under the item labeled "z ~ 1 composite" ) . 

Our z < 0.05 sample of early-type cluster galaxies, 
selected in the same manner as our z > 0.3 sample, has 
a SO fraction of 68 ± 3%. This fraction is consistent with 
other analyses. Examining the whole sample of early- 
type galaxies within 2R2oq/'k, we find a SO fraction of 
67±3% for all of the early- typ e galax ies within for the 10 
clusters we use from Dressier (1980a), in good agreement 
with that found by Dressier et al. (1997) . Therefore, our 
red-sequence selection produces a comparable sample of 
early-type cluster galaxies to those from previous efforts. 

We find no statistically significant evidence for evolu- 
tion within the z > 0.4 cluster sample, whi ch is also con- 
sistent with most of the previous work (e.g.. | Fasano et al.l 
l2000HPostman et al.ll2005HDesai et al.ll2007[ ). Using the 
visual classifications, our sample also shows that evolu- 
tion in the SO fraction occurs between z ~ and z ~ 0.4 
as the authors above have noted. Note that the fractions 



we are discussing here are the fraction of SO galaxies 
in th e early-type galaxy population (see iPostman et al.l 
l2005l for a discussion of the size of the systematic error 
in separating the elliptical and SO populations.) 

The result from our measurements of the ellipticity dis- 
tributions in our clusters is that there is essentially no 
evolution in the ellipticity distributions from z ~ 1 to 
z ~ 0, and thus that there is no change in the over- 
all distribution of the bulge-to-disk ratio of early-type 
galaxies. If we assume no evolution in the bulge-to- 
disk ratio distributions of the elliptical and SO popula- 
tion separately as others have done, we can conclude 
that the E/SO ratio does not evolve over this inter- 
val. This contrasts with the clear and significant evo- 
lution Jn_the SO fractio n seen from ma n y studies (e.g.. 
Dress ier et al.l [1 997: Fa sano et alj l2000t [Postman etall 
120051 : iDesai et alJ 120071) based on visual classifications. 
While one explanation is that there is a problem with 
th e visual classific a tions, this raises an int e restin g issue. 
In IPostman et al.l (|2005D and iDesai et al.l (|2007D . what 
evolves at z > 0.4 is the fraction of SO galaxies and 
the fraction of spiral and irregular galaxies. In contrast, 
the fraction of elliptical galaxies does not evolve. If the 
misclassification of round SO galaxies as ellipticals alone 
causes the fraction of SO galaxies to decrease, it is puz- 
zling that there is no corresp onding increase i n the frac- 
tion of elliptical galaxies in IPostman et al.l ()2005l ) and 
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Fig. 5. — Median ellipticity versus redshift for the clusters in our 
sample. The open circles are the median ellipticities for the early- 
type galaxies in each cluster, after applying our magnitude and ra- 
dius selection. The median ellipticity, e^^j = 0.29 ±0.02, for all of 
the galaxies in z < 0.05 cluster sample is shown by the blue square 
and line. The other squares are the medians for all of the cluster 
galaxies in the redshift ranges 0.3 < z < 0.6 (green), 0.6 < z < 1.0 
(orange), and z > 1 (red) samples, respectively. The median ellip- 
ticity for the whole of the z > 0.3 sample is e^ed = 0.30it0.01. The 
median ellipticity for cluster galaxies in the range of 0.3 < z < 0.6 
is 0.31 it 0.02, while the median ellipticity for cluster galaxies in 
the range of 0.6 < 2 < 1.0 is 0.30 ± 0.02. At z > 1, the me- 
dian ellipticity is Emed = 0.29 ±0.03. All values are all in excellent 
agreement with the 0.29 ±0.02 we find for the z < 0.05 sample. We 
find no individual clusters that have drastically different ellipticity 
distributions. The lack of any trend in the ellipticity, and minimal 
cluster-to-cluster variance, is a striking result. The median ellip- 
ticity at z > 0.3 is statistically identical to that at z < 0.05, being 
higher by only 0.01 ± 0.02. 



iDesai et al.l (|2007D . Rather, in previous work, the spi- 
ral fraction increases, and there is httle expectation that 
face-on SO galaxies would be misclassified as spirals. The 
source of this difference remains to be resolved. Our 
study does not provide an opportunity for a resolution 
of this issue, but we wanted to highlight it for further 
work by others. 

5.2. Ellipticity by Galaxy Type 

Our simulations show that the estimates we use for el- 
lipticities are robust, and do not show a large systematic 
change with redshift (see Figure[TJ) Some of the observed 
change in the fraction of SO galaxies in Figure [7] may 
come about from a mis classification of round ga laxies 
as ellipticals, as others pOrgensen fc FraiD3ll994D have 
suggested and as we have discussed earlier. It would be 
useful to examine a sample of SOs for which the elliptical 
contamination was, likely, very small. We examined the 
fraction of galaxies with ellipticities above the median 
of the z < 0.05 SO population, €,ned = 0.38. Measur- 
ing the fraction of galaxies with emed > 0.38 should give 
us an estimate of the SO population that we expect to 
be modestly contaminated by elliptical galaxies. In our 
z < 0.05 sample, we find that 5 ± 3% of galaxies with 
e > 0.38 are ellipticals. We plot the fraction of galaxies 
with emed > 0.38 in Figure [5] as a function of redshift. 
This once again, shows no change with redshift, imply- 
ing a SO fraction that does not evolve under the assump- 
tion that the ellipticity distribution of ellipticals and SO 
galaxies do not separately evolve. 



If the SO fraction decreases with redshift while the 
median ellipticity of the early-type population stays the 
same, we expect that the ellipticity distribution for the 
SO galaxies should increase. In Figure [H we plot the me- 
dian ellipticity of galaxies that are classified visually as 
SO galaxies. We show that the median ellipticity of the 
SO population is higher in the higher redshift clusters. In 
the z > 0.3 sample, we find that the median ellipticity of 
the SO galaxy population is emed = 0.47 ± 0.02 as com- 
pared with Emed = 0.38±0.02 in our z < 0.05 low redshift 
sample^^. The difference between the z > 0.3 and the 
z < 0.05 SO samples medians is significant at the > 3cr 
when using a Student's t-test. Restricting the redshift 
range to z > 0.4, the redshift where other authors have 
reported that the SO fraction significant evolves, we find 
emed — 0.48 ± 0.02, which is also a > 3(t difference. The 
significance remains > 3a regardless of whether Coma is 
included in the z < 0.05 sample or not. 

In Figure \W\ (similar to Figure [H but with just the 
complete sample with E and SO identified), we plot the 
ellipticities as a function of absolute magnitude in four 
redshift bins. It appears from this figure that there is 
a deficit of rounder SO galaxies (blue points) at higher 
redshifts (excluding the z > 1 bin which is less complete 
and has poor statistics), implying that the SO galaxies 
at z > 0.3 are drawn from a different distribution than 
those at z < 0.05. The progression of higher ellipticities 
for the SO population thus appears to come about from 
a change in the shape of the ellipticity distribution of 
SO galaxies at higher redshift as compared with those at 
lower redshifts. We discuss this in H5.4\ 

In Figure we quantify what we see in Figure 1101 
We plot the distributions of the ellipticities for the z > 
0.3 elliptical and SO galaxies separately. We find that 
the shape of the z > 0.3 SO distribution differs from 
the z < 0.05 distribution, in contrast with the elliptical 
distributions which show little change between the two 
samples. For example, if we examine the number of SO 
galaxies with e < 0.3, we find that there are 56 in the 
z > 0.3 sample, whereas we expect 95 from the scaled low 
redshift sample of SO galaxies. We compared the shapes 
of the two SO ellipticity distributions using a number of 
statistical tests. The probability of such a change in the 
ellipticity distribution of SO galaxies is > 3(T when using 
a Kolmogorov-Smirnov, a Wilcoxon-Mann- Whitney rank 
sum test, and a Kuiper test. These are the same tests 
that showed no evolution in the ellipticity distribution in 
§4.21 These tests reinforce the result, from the analysis 
of the trends in Figure [9l that the cluster SO galaxies 
have a different ellipticity distribution at z > 0.3 than at 
z < 0.05. These trends, both in the median and in the 
shape of the distribution, are present regardless of how 
we select the higher redshift sample, whether z > 0.3, 
z > 0.4 or 0.4 < z < 1. In all cases, the distributions 
differ from the z < 0.05 sample at > 3ct. 

Our work is not the only one that finds a different me- 
dian elliptici ty for SO galaxies at higher redshift in clus- 
ter samples. iMoran et al.l (|2007f l give the morphological 

These median values are the medians of the all SO galaxies 
that meet our selection criteria. These are not the median of the 
data points we plot in Figure [9] The error we quote is the error on 
the median, or ^J\it/2) the error on the mean. We confirm these 
errors with bootstrapping which yields good agreement. 
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Fig. 6. — Cumulative (top) and differential (bottom) ellipticity distributions for our four rcdshift bins. In each plot, we show the 
distribution of the ellipticities of all early-type galaxies in our magnitude and radius selection for all of the cluster galaxies in a given 
redshift bin. The distributions are represented by the shaded histograms, outlined in blue for 2 < 0.05, green for 0.3 < z < 0.6, orange for 
0.6 < z < I, and red for z > 1. For comparison with the z > 0.3 cluster samples, we plot our sample of z < 0.05 early-type galaxies in blue, 
with the low-redshift line "hidden" when the two lines (frequently) overlap. The z < 0.05 sample is normalized to have the same number 
of galaxies as each of the z > 0.3 cluster samples. The small deficit of highly elliptical galaxies at z > 1 is not statistically significant 
because of the small sample size. Again, it is striking that there is no evolution in the shape of the distribution from 2 < 0.05 to 2 > 1. 
The distribution of ellipticities at z > 0.3 agrees with the shape of the z < 0.05 distribution at the 1-2% level (i.e., the probability that they 
are drawn from the same distribution is 98-99%). 




Fig. 7. — Ratio of m orpho logically identified SO galaxies from 
the studies discussed in f\2.2\ to the total number of elliptical and 
SO galaxies. The fraction of SO cluster galaxies at 2 > 0.4 is lower 
than that seen at 2 < 0.05, as seen in previous work. We compute 
the average fraction of SO galaxies in the early-type population for 
all of the 2 > 0.3 clusters and plot that value, 42 ± 2% as a solid 
black line. For contrast, we show the 2 < 0.05 value, 68 ± 3% , as 
a solid blue line. We plot, as squares, the average values of the SO 
fraction of early-type galaxies in three redshift bins, 0.3 < 2 < 0.6 
(green), 0.6 < 2 < 1.0 (orange) and 2 > 1 (red). Our sample shows 
the same trend in SO fraction found by other work. However, the 
very different fraction of galaxies classified as SO's at 2 > 0.4 is in 
seeming contradiction with the lack of evolution in the ellipticity 
distribution. 

classifications and ellipticities of a sample of SO galaxies 
at z ~ 0.45. The authors measured the ellipticities of 
the cluster galaxies using the same approach and soft- 



FlG. 8. — Fraction of cluster early-type galaxies with ellipticities 
greater than the median ellipticity of the 2 < 0.05 SO galaxies. 
The open circles are the fractions of galaxies with tmed > 0.38, 
in each cluster. The squares are the fractions for the 2 < 0.05 
(blue), 0.3 < 2 < 0.6 (green), 0.6 < 2 < 1.0 (orange), and 2 > 1 
(red) samples respectively. At low redshift, this fraction has only a 
small contamination from elliptical galaxies, such that the fraction 
of galaxies with Emed > 0.38 is the half fraction of SO galaxies with 
a 5% contamination from elliptical galaxies. If the observed lack 
of evolution in the median ellipticity of the 2 > 0.3 population is a 
result of evolution in the ellipticities of elliptical galaxies masking 
the decline in the SO population, we should still see a change in the 
population of galaxies with eme^i > 0.38. Galaxies with e^^j > 
0.38 represent the most inclined of the disk-dominated early-type 
galaxies at low redshift. The lack of evolution in this fraction 
suggests little evolution in the fraction of disk-dominated early- 
type galaxies as a whole. 
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Fig. 9. — Median ellipticity versus redshift for cluster SO galaxies. 
The open circles are the median ellipticities for the galaxies in each 
cluster. The blue squares are the median for all of the galaxies in 
the z < 0.05 cluster sample. The other squares are the medians 
for the 0.3 < z < 0.6 (green), 0.6 < z < 1.0 (orange), and z > 1 
(red) samples, respectively. We also show the median value for 
the z < 0.05 SO sample as a blue line. We show, with a black 
line, the median ellipticities of all the SO galaxies in our z > 0.3 
sample. The median ellipticity of the SO population increases at 
z > 0.3, with emed = 0.47 ± 0.02. The median ellipticity of the 
z < 0.05 SO galaxy is Cmed = 0.38 ± 0.02. Such a change in the 
ellipticity of galaxies classified as SOs can explain how the fraction 
of galaxies classified as SO's decreases while, at the same time, the 
median ellipticity of the cluster galaxy population as a whole stays 
the same. 

ware that we use (and thus their eUipticit ies are PSF- 
corrected.) For their cluster SO galaxies, iMoran et al.l 
([^07) find emed = 0.48 ± 0.04, in excellent agreement 
with our value of emed = 0.48 ± 0.02. 

Previous work found little or no evolution in the ellip- 
ticities of SO and elliptical galaxies. However, most of 
the previous work relied on measuring the ellipticities of 
galaxies using flux- weighted moments of galaxies, with- 
out taking into account the "blurring" effects of the PSF 
or the lower resolution of observations at higher redshifts. 
We show, in Appendix [Al that the ellipticity measure- 
ments of galaxies based on simple flux- weighted moments 
are biased quite significantly towards rounder values for 
the ACS and WFPC2 data that has been used without 
correcting for the effects of smoothing by the PSF. 

5.3. Constraints on the Evolution of the Bulge-to-Disk 

Ratio 

The observed lack of any evolution in the ellipticity 
distribution for early-type galaxies has interesting impli- 
cations for the distribution of the bulge-to-disk ratios of 
the early- type population from 2; ~ to z ^ 1. The aver- 
age observed ellipti city is directly related to t he average 
intrinsic ellipticity (jBinnev fc Merrifieldl[l998h . If the av- 
erage bulge-to-disk of a population of galaxies changes, 
the average intrinsic ellipticity must change too, and vice 
versa. Our measurement of an unchanging ellipticity dis- 
tribution implies that the average bulge-to-disk remains 
constant. 



We performed simple simulations to derive the aver- 
age ellipticities using models of galaxies with a variety 
of bulge-to-disk ratios projected over all possible view- 
ing angles with a fixed scale length for the bulges and 
disk. We find a roughly linear relation between the av- 
erage bulge-to-total ratio and the average ellipticity of 
the population. For example, for these simple models, 
a change in the average bulge-to-total ratio from 0.1 to 
0.3, or from 0.3 to 0.6, would be observed as a shift of 
5e — 0.08 in the average ellipticity. Such a shift would 
be a 4(7 change for our samples (see Figure [5l) While 
one could potentially develop a model whereby the un- 
derlying scale lengths relatively evolve with redshift, this 
seems rather contrived. "Occam's Razor" would lead 
one to prefer a model where a constant ellipticity dis- 
tribution implies a constant distribution of bulge-to-disk 
ratios over the redshift range from z ^ 1 to z ^ 0. 

5.4. Implications for Evolution in the Fraction of 
Disk- Dominated Early-type Galaxies 

We performed simulations to see if the observed SO 
fraction of 42%, for the sample of clusters with z > 0.4, 
is consistent with the lack of evolution we see in the el- 
lipticity distribution. We constructed 10,000 mock cat- 
alogs for the clusters in the z > 0.4 sample. For each 
cluster, we assumed that the fraction of SO galaxies was 
42%. We randomly drew the ellipticities for the SOs from 
the observed distribution of SOs at z < 0.05, and drew 
the ellipticities for the Es from the observed low red- 
shift distribution for Es (see Figure [3l) We added the 
uncertainties and systematics applicable at high redshift 
as detailed in §3.11 and Appendix IA.2I The resulting el- 
lipticity distributions had median ellipticities lower than 
the observed value (0.29) for all 10,000 of the mock cat- 
alogs, with a typical value of emed = 0.25. This rules out 
at > 4cr the assumption that the underlying ellipticities 
and SO fraction are both constant with redshift. 

Our others tests which we use to compare the z < 0.05 
sample and the z > 0.3 sample give similar results 
(the Kolmogorov-Smirnov statistic, a Wilcoxon-Mann- 
Whitney rank sum statistic and a Kuiper statistic). None 
of the mock catalogs show the same good agreement be- 
tween the z > 0.3 sample and z < 0.05 sample in the 
ellipticity distribution as we observe. 

We repeated these simulations for a variety of SO frac- 
tions. We find that the SO fraction value that best 
matches the z > 0.4 ellipticity distribution is 66% with 
a 1(7, or 68% confidence limits, of ±6.5%. Another way 
to consider this is a 1(7 change in the median ellipticity 
at z > 0.3, a change of 5e = 0.02, would imply a change 
in the SO fraction of 6.5%. We can rule out, at the 95% 
confidence limit, SO fractions below 53%, and we can rule 
out the observed fraction of visually classified SO galaxies 
of 42 ± 2% at the > 99.9% confidence fimit or at the 4cr 
level. 

To reconcile the evolution seen in the SO fraction in the 
visually-classified samples with the lack of evolution in 
the ellipticity distribution requires that either the SO and 
elliptical populations evolve in ellipticity in such a way 
that the combined samples shows no evolution, or that 
some fraction of the SO poptflation has been misclassi- 
fied as other morphological types (cf.. .j0re.'cnscn fc Franj3 
I199I iBlakeslee et al.ll2006D . 
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Fig. 10. — The distribution of ellipticity versus Mg for cluster elliptical (red) and SO (blue) galaxies in four redshift bins. The solid 
points are our sample within our magnitude limits of —19.3 > Mg + 1.208^ > —21. This plot is similar to Fig. |4] except that we color 
code the SO and elliptical galaxies separately. We show the median values of the populations as lines, red for the ellipticals and blue for 
the SO galaxies. Wc find no statistically significant evolution in the elliptical population. There is a statistically significant (> 3(t) change 
in the distribution of ellipticities for SO galaxies when comparing the z < 0.05 SO population with Cmed = 0.38 ±0.02 and the 0.3 < z < 0.6 
{^med = 0.46 ± 0.02) or 0.6 < z < I sample (tmed = 0.49 it 0.02). The uncertainty in the highest redshift sample is too large for any useful 
comparisons to be made. The change in the ellipticity of the elliptical galaxies with redshift is not statistically significant. From this figure 
we conclude that the increase in the median ellipticity seen in Fig. [9] could arise from a lack of round SO galaxies at redshifts higher than 
z > 0.3. 
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Fig. 11. — Ellipticity distributions for cluster elliptical (left) and SO (right) galaxies. The z < 0.05 sample is shown in blue, while the 
z > 0.3 sample is shown in red, with the median ellipticity for each shown as a appropriately colored vertical line. We find that the z > 0.3 
elliptical population is statistically the same as the z < 0.05 population. In contrast, there is a statistically significant difference (> 3(t) in 
the ellipticities of the SO population. This appears to be a deficit of round SO galaxies. We find 56 SO galaxies with e < 0.3 in our z > 0.3 
sample, whereas we expect 95 SO galaxies, after rescaling the z < 0.05 to match the total number of expected SO galaxies. This reinforces 
Fig. O showing that the typical ellipticity for SO galaxies at 2 > 0.3 is higher than at 2 < 0.05 and it appears that there is a deficit of 
round SO galaxies which is causing this difference in the ellipticity distributions. 
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5.4.1. Possible Spiral Contamination? 

One possibility is that the higher eUipticity, or more 
edge-on, SO population is contaminated with misclassi- 
fied red spirals. This would cause the SO fraction to be 
overestimated, while, at the same time change the shape 
of the eUipticity distribution. The spiral fraction among 
galaxies in the ma ss range of the clus ter sample in our 
work is ~ 10-15% (jHolden et al.ll2007D while the f raction 
of dusty, red objects is ^ 10% (jSaintonge et al.l [20081 ) . 
Since the misclassification would be predominately for 
edge-on systems, these fractions are likely upper limits 
for the fraction of spirals misclassified as SO galaxies. Is 
it possible that such a small fraction of "spiral" galax- 
ies could contaminate the population visually classified 
as SO galaxies and change the eUipticity distribution at 
higher redshift? We examined the ellipticities of the 
galaxies classified as spirals in both the z < 0.05 and 
z > 0.3 sample to see if there is a deficit of edgc-on 
spirals. We found the opposite, that there is a shift in 
the high redshift spiral population to larger ellip ticities, 
Cmed = 0.38 i 0.04 at z < 0.05 (as expected, see iRvdenI 
120061 ) to e^ed = 0.45 ± 0.03 at z > 0.3, but it is not sta- 
tistically significant. From this we conclude that a red 
spiral population is not a major source of contamination, 
especially at large ellipticities. 

Given the quantitative nature of the current study we 
feel that it provides a key datum for consideration of 
the evolutionary history of early-type galaxies over the 
last 8 Gyr. The result that the eUipticity distribution 
is essentially unchanged from z ~ to z ~ 1, combined 
with the result that the fraction of highly elongated disk- 
dominated systems does not evolve (see Figure El) , pro- 
vides strong evidence that the overall bulge-to-disk ratio 
distribution of the population does not evolve. 

6. SUMMARY AND IMPLICATIONS 

We have compiled a sample of 10 z > 0.05 clusters of 
galaxies and a comparable sample of 17 z > 0.3 clus- 
ters of galaxies with HST imaging. For each cluster, we 
selected a subsample of galaxies that lie on the red se- 
quence and have been classified as an early-type galaxy 
(elliptical or SO galaxy). To derive a robust sample at 
all redshifts, we selected galaxies to fall within a magni- 
tude range of —19.3 > Mb + 1.208z > —21 and required 
the galaxies to lie within 2i?2oo/7'' (to ensure that the 
de-projected sample lies within i?20o) of the cluster cen- 
ter. We change the rest-frame B magnitudes by 1.208z 
to match the mass-to-light evolution of cluster early- 
type galaxies as measured by the fundamental plane 
(Ivan Dokkum fc van der Marell[200l . The magnitude 
limits correspond to M^ + l > Mb > M^-0.75 or galax- 
ies with stellar masses roughly between 10^°'^ Mq < 
M < 10"-2 Mq, assuming a "diet" Salpeter IMF. 

We performed extensive simulations of the eUipticity 
measurements to test their robustness and to minimize 
systematics with redshift. We found that it is crucial 
to use model fits (e.g., GALFIT) where the effect of the 
PSF can be modeled and removed. From our extensive 
simulations we find that we can robustly measure the 
eUipticity of cluster early-type galaxies (by using PSF- 
corrected model fits) out to z ~ 1.3 at magnitudes that 
allow the galaxies to also be morphologically classified. 
The systematic error we find in our simulations is only 



(5e = —0.01, ora3% change for the median galaxy. We do 
not apply this small correction to our data. We find that 
not including the effect of the PSF can cause systematic 
shifts of 6^ — —0.10, or a 30% change for the median 
galaxy. 

Using our two samples of cluster early- type galaxies, 
we measure the evolution in the eUipticity of z > 0.3 
cluster early-type galaxies compared with a substantial 
sample at z < 0.05. 

1. We find no evolution in the median eUipticity of z > 
0.3 cluster early-type (eUiptical and SO) galaxies 
relative to the low redshift sample. The median 
eUipticity at z < 0.05 is €med = 0.29 ± 0.02, and 
is emed = 0.30 ± 0.01 at z > 0.3. The median 
eUipticity of our sample of J^Sl early-type galaxies 
at z > 0.3 is statistically the same as that of the 
210 early-type galaxies at z < 0.05, being higher by 
only 0.01 ± 0.02 or 3 ± 6%. 

2. The shape of the eUipticity distribution of z > 0.3 
galaxies also does not evolve. The distribution of 
ellipticities at z > 0.3 agrees with the shape of the 
z < 0.05 distribution at the 1-2% level (i.e., the 
probability that they are drawn from the same dis- 
tribution is 98-99%). 

3. Using visual classifications from previous reference 
studies, we find that our sample shows a similar 
decrease in the fraction of early-type galaxies clas- 
sified as disk-dominated systems (SOs) at z > 0.4. 
As other studies have found, t he SO fraction is 
68±3% at z < 0.05 and decreases to 42±2% at 
z > 0.4. 

4. For the SO fraction to decrease with increasing red- 
shift, while the median eUipticity of the overall clus- 
ter early-type population stays the same, the me- 
dian SO eUipticity needs to be larger at higher red- 
shift. There is a trend in our data for z > 0.3 in 
this sense. The median SO galaxy at z < 0.05 is 
Emerf = 0.38 ±0.02 while at z > 0.3, we find a. ^ ia 
different emed — 0.47 ± 0.02 for galaxies classified 
as SOs. Whether this change is large enough to ac- 
count for the change seen in the morphological sam- 
ple remains to be determined. We explored if the 
change in the eUipticity distribution of SO galax- 
ies could come from a deficit of round SO galaxies 
or from the addition of misclassified spiral galax- 
ies. However, we found that neither of these simple 
scenarios can explain our results. 

Our results on the unchanging eUipticity distributions 
lead us to conclude that there has been little or no evo- 
lution in the overall distribution of bulge-to-disk ratio of 
early-type galaxies over the redshift range < z < 1 
for morphologically-selected samples of early-type red- 
sequence galaxies with -f 1 > Mb > Mg — 0.75 in 
the dense cores of clusters inside i?200- In particular, 
our results allow us to rule out a SO fraction of <47-51% 
at z > 0.3 at the 3cr level assuming no evolution in the 
eUipticity distributions of elliptical and SO galaxies. // 
we assume, as in all previous studies, that the intrinsic 
eUipticity distribution of both elliptical and SO galaxies 
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remains constant, then we conclude from the lack of evo- 
lution in the observed early-type ellipticity distribution 
that the relative fractions of ellipticals and SOs do not 
evolve over the last ^ &Gyr, or from z ^ 1 to z — 0, 
for a red-sequence selected sample of early-type galaxies 
in the cores of clusters of galaxies. 

These results do highlight an inconsistency within the 
wide range of studies that have occurred on the evolution 
of early- type galaxies over the last decade. Over this red- 
shift range, and particularly since z ~ 0.4, the fraction of 
morphologically-identified SO galaxies from visual classi- 
fications has been foun d to be significantly l ower than the 
z ^ value (see, e.g ., IPressler e t al.lll997l: iFasano et al.l 
120001 : iPostman eta l. 2005; Dcsa i et al.l 120071) . Reconcil- 
ing the trends seen in the elliptical, SO and spiral frac- 
tions, as discussed in !j5. II remains to be understood. 

The early-type galaxy population, both in the field 
and in clusters, does not evolve purely passively. The 
volume-averaged number densi ty of red galaxies , which 
are mostly early-type galaxies ()Bell et alll2004aD, grows 



by a factor ~ 2 between -z ^ 1 and the present (|Bell et al.l 
I2004U iBrown et all l2007t iFaber et al.l 120071) . For the 
cluster popul ation such evolution is harder to quantify , 
but mergers (|van Dokkum et al.|[l999l : iTran et al.|[2005f) 

and filaments aroimd massive clusters at z ~ 0.8 (e.g., 
iKodama et al. 2005; Patcl et al. 200i) suggest that in- 
teractions and infalling galaxies enhance and modify over 
time the galaxy population in the cluster core. Besides 
evolution in the population of early- type galaxies, in- 
dividual early-type galaxies also change over time. As 
pointed out by, e.g., |j0r gensen et al. (2005), the evolu- 
tion of the line strengths is not compatible with purely 
passive evolution. Moreover, early-type galaxies at z ~ 1 
were recently demonstrated t o be significantly smalle r 
than today by a factor of two Ivan der Wei et aU (|2008D . 
Hence, the resulting picture is complicated: early-type 
galaxies, both as individual objects and as a popula- 
tion, undergo substantial changes between z ~ 1 and 
the present. 

On the other hand, many basic properties of the early- 
type population have remained unchanged over the past 
~ 8Gyrs, providing useful means for characterizing and 
quantifying the observed evolution. At any redshift 



z ^ 1 early-type galaxie s occup y a tight color-magnitude 
relation (e.g.. Blakcslce et al.l [2006) and fundamental 
plane (e.g.. IWuvts et al.ll2004f ). which suggests smooth 
and regular evolution. Moreover, neither in the field 
nor in clusters has the early-type fractio n changed sig- 
nificantly in mass-selec ted samples ( H olden et al.l [20071 : 
Ivan der Wei et aLll2007f ) . Together with these character- 
istics, the results presented in this paper fit into a picture 
in which the field and cluster early-type galaxy popu- 
lations grow and change, but only while leaving many 
basic characteristics the same. In particular, we have 
demonstrated from the constancy of the ellipticities of 
early-type galaxies over the last ~ 8Gyr that the bulge- 
to-disk ratio distribution of the cluster population re- 
mains constant. This suggests that processes that change 
the bulge-to-disk ratio of individual early-type galaxies 
and the bulge-to-disk ratios of newly formed or accreted 
early-type galaxies are balanced such that the overall 
bulge-to-disk ratio distribution remains the same. It re- 
mains an open question why this is the case, but it is 
clear that much can be learned about the formation pro- 
cess of early-type galaxies by studying their properties 
in even more detail and extending their observation to 
higher redshifts. 
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APPENDIX 

A. ROBUSTNESS OF ELLIPTICITY MEASUREMENTS 

We made a number of additional tests of the robustness of the ellipticity measurements with the goal of estimating 
the magnitude of any systematic errors. These include adding fake galaxies into our imaging data, comparing the 
ellipticity measurements made using SExtractor to those using GALFIT, and using the a range of (incorrect) PSFs to 
estimate the possible sensitivity to systematic errors. 



Al. Simulations with Fake Galaxies 

We performed a number of simulations using fake galaxies to test the reliability of the ellipticity measurements. 
These simulations were done in order to establish the best approach to use and to assess how well we can measure the 
ellipticity, given our data. 

The artificial galaxies were made by simulating two Scrsic models with the same centroid, one with rt = 4 and one 
with n = 1. The two components had the same luminosity, but different effective radii. The first effective radius was 
chosen to be three times the resolution of the ACS camera, while the second was nine times the ACS resolution. Each 
component was given the same ellipticity. These represent well-resolved lenticular galaxies. These fake images were 
then simulated at magnitudes typical of the early-type galaxies in the z 0.8 — 0.9 clusters in our ACS imaging, and 
then added to a real ACS image. Reassuringly, the simulations recovered the input ellipticities to a high accuracy, 
with a scatter of just CTg/e = 0.027 and an offset of Se/e ~ —0.0007. The scatter and the offset was highest for the 
most round objects. For objects with an input ellipticity of 0, the average recovered ellipticity was 0.03 with a scatter 
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of CTe = 0.014. As expected, objects that would be seen as perfectly round given an arbitrarily high signal-to-noise, 
are not measured to be round at typical signal-to-noise levels, but the offset, as noted above, was quite small. 

These simulations were made for a variety of bulge-to-disk ratios, by varying the relative luminosities of the two 
components. In addition, we assumed that the intrinsic ellipticity of the n = 4 bulge was e = 0.3, or e = 0.2 when 
averaged over all projections or over all cos(i). 

A2. Simulations with Real Galaxies 

We compiled a number of images of low redshift elliptical and SO galaxies. We began with the catalog of iFrei et al.l 
(|1996( l. That catalog contains 24 elliptical and SO galaxies, including two classified as SO/a. We expanded upon this 
by including 49 additional ear ly- type galaxies in the Virgo cluster. Each was selected to be in the Virgo cluster sample 
discussed in iMei et al.l (|2007f ). Most of the galaxies in that sample have accurate distances from surface brightness 
fluctuations. Those that do not have distances were assumed to be at the center of t he main V i rgo co ncentration. 

Once we selected the list of early- type galaxies, we extracted images of those not in lFrei et al.l (|1996f l from the SDSS 
DR5 in the g filter. For those images from the SDSS, we removed nearby bright stars and other galaxies, replacing 
the pixels with those randomly selected from nearby blank regions. This produced a final collection of 73 images of 
nearby galaxies imaged in either g ov B. 

We added these galaxies, appropriately scaled in size and magnitude, to the imaging data used for each of the 
clusters, and replicated the detection and measuring process. The galaxies were resized as appropriate for the redshift 
of the cluster, convolved with a representative PSF, and then noise was added to each. Each galaxy was simulated at 
a variety of input magnitudes, and multiple realizations w e re do ne at each magnitude to evaluate the effect of noise 
and binning. This is the same process used in lHolden et all (|2004f) . The magnitudes were selected to span the range of 
observed galaxy magnitudes, to 0.4 mag below the limit used for classifying galaxies in each sample. This is important 
to ensure that the selection at the limit is not influenced by the magnitude cutoff. 

Recovery of the Input Galaxy's Ellipticity 

We compared the offset between the ellipticity measurement from the original image and that of the simulated high 
redshift image. In general, the offset between the low redshift original image and high redshift scaled image ellipticities 
are consistent. For galaxies within our magnitude limits, we found that the typical systematic offset is from (5^ —0.03 
to (5e ~ 0.01, where the simulated galaxies are measured to be less elliptical at high redshift than in their original, low 
redshift images. For galaxies at or below our selection limits 5^ ^ —0.06 to 5^ ~ 0.05. The median systematic offset 
was 5^ ^ 0.01 for our whole sample of simulated galaxy images. We show the ellipticities measured in the output of 
the simulations as a function of the input ellipticities in Figure [TJ The scatter scales with the input signal-to- noise as 
expected. One exception is that, in general, the measurements made with WFPC2 are worse than those made with 
ACS at the same redshift. This is expected, because WFPC2 is under-sampled and the images are typically at a lower 
S/N. 

Random Errors 

We used the multiple realizations of the galaxy simulations to examine the random errors. The scatter in the 
measurements for a given input galaxy in a given cluster observation is, on the whole, small. The ACS imaging has 
measurements with a median scatter of ~ 0.01 in ellipticity, while the WFPC2 imaging shows larger median scatter 
of (Te ~ 0.02 — 0.03 in ellipticity. These numbers are close to the typical errors reported by GALFIT. GALFIT's reported 
errors scale with the magnitude of the galaxy, which we also find in our simulations. At the faintest magnitudes, at 
or below the magnitude limits for classification, we find a scatter on the ellipticities as high as ^ 0.05, once again 
matching the scatter in the ellipticity as estimated by GALFIT. 

Systematic Errors from PSF Mismatches 

Because the PSF so important, we performed simulations where GALFIT was run using a different PSF than the 
PSF that was used to convolve with the data. We reran the simulations for MS 1054-03 using the correct 1775 PSF 
when simulating the z — 0.83 galaxy, b ut we ran GALFIT u sing the /814 PSF. The differences between all three PSFs 
in ter ms of encircled energy are small (ISirianni et al.l 120051 ). However, there is an asymmetry in the PSF of the 
filter (|Sirianni et al.|[2005l : iJee et al.ll2007[ ) that causes a noticeable difference in the resulting PSFs. 

We show the results of using the /814 PSF to fit in Figure [12] We find an overall offset of 5^ ^ 0.027 between the 
average ellipticities, which is of order the same size as the bins we use in the top panel. Examining the lower panel, we 
can see that the offset is not constant with ellipticity, but also depends on the input galaxy used for the simulation. 
For example, along the x-axis at e ~ 0.3, there are a number of points where the incorrect PSF yields e ~ 0.4 while 
other galaxies have similar ellipticities when using the correct or incorrect PSF. From this we can conclude that typical 
systematic errors on the ellipticity measurements are on the order of 5^ ~ 0.01 — 0.03, but this depends on the shape 
of the galaxy. The PSF is more asymmetric in the redder passbands, so the systematic error is likely to be largest 
for the highest redshift clusters. The statistic we use the most is the median ellipticity. Using these simulations, we 
find that the offset in the median ellipticity is 5^ = 0.018 for the Isi4 data. This is a quite a small effect, given the 
noticeable changes that we have made to the PSF. 
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Fig. 12. — EUipticity determined with wrong PSF versus ellipticity with the correct PSF. The bottom panel shows the output ellipticity 
as a function of the input eUipticity, with a line with a slope of one to guide the eye. The top panel shows the input ellipticity distribution 
in blue and the output as the gray-shaded histogram. The simulations were done for cluster galaxies at z = 0.83 which were convolved 
with the 4775 PSF, but the modeling of the galaxy luminosity profiles was done using the /814 PSF has some asymmetry. The use of the 
incorrect PSF in this case has only a modest impact, biasing the ellipticities towards higher values by 5^ ~ 0.03 at the highest ellipticities 
and becomes smaller, 5^ = 0.018 or a ~6% change, at the median input ellipticity. 

A3. Comparison of Ellipticity Measures 

Previous work, such as lDressler et al.1 (|1997f l and iPostman et all (|2005f ). use SExtractor's flux- weighted moments to 
estimate the ellipticities of galaxies. This approach, however, does not remove the effect of the PSF. Below we will 
compare the ellipticity measurements we use with those used by previous work in the literature. We will also use our 
simulations to evaluate how well these techniques reproduce the underlying ellipticity distribution. 
In FigurefO} we compare our ellipticity measurements for CL 0016-1-16 CL 0024-1-16 and CL 0054-27 with those from 



iSmail et all ([l997h . Our measurements, as we stated above, are based on fitting e lliptical models of Sersic profiles 
to the data, including t he effects of the PSF. The results from ISmail et"aI1 (|1997t l are based on using SExtractor 
(jBertin fc Arnoutill996l ). SExtractor estimates the ellipse by computing the flux- weighted second-order moments of 
the gal axy. SExtra c tor in cludes those pixels above the threshold defined by an isophote as measured in a smoothed 
image. ISmail et"al] (|1997f) used a detection isophote of ~ 1.3cr and a O'.'S diameter, or 3 pixels with WFPC2, top-hat 
filter for the smoothing. 
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Fig. 13. — EUipticity measurement comparison with ISmail et ahl l|1997). Our ellip ticities measured by fitting PSF convolved models to 
the images with GALFIT are on the x-axis. The y-axis values are ellipticities from ISmail et al.l II1997I ) measured from the second-order 
fiux-weighted moments of galaxies using SExtractor. Galaxies classified by eye as SOs are open blue squares while ellipticals are red circles. 
The smaller symbols are galaxies with radii less than 0"3, or ~ 2 WFP C2 resolution eleme nts. The black lino is a line of slope one with an 
intercept of zero. Our ellipticities are generally less round than those of lSmail et al.l (11997'), and this is discrepancy is larger for the smaller 
galaxies. The median offset in the ellipticities for SO galaxies is 5trned = —0.13, much larger than the worst systematic error we found in 
our simulations of our ellipticity measurements. This highlights the importance of including the PSF in the measuring process. 



We find that our ellipticities are less round than those from lSmail et al.l ()1997l ). We hypothesize that this comes about 
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for two reasons. First, our approach removes the smoothing from the PSF which tends to circularize the eUipticity 
measurements. Second, much of our data was observed with ACS, which has a more compact PSF than WFPC2. The 
median offset in the eUipticities for SO galaxies is Scmed = —0.13. This is a very large change, much larger than any 
of the systematic error we found in our simulations of our eUipticity measurements. This highlights the importance of 
including the PSF in the measuring process, and does indicate the challenge of using WFPC2 data for measurements 
at higher redshift. 
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Fig. 14. — SExtractor measures of the eUipticity versus those from GALFIT for ACS imaging data. This plot is similar to Fig. 1131 and 
uses the same symbols, but compares eUipticities from measurements of the second-order flux-weighted moments made with ACS imaging, 
as opposed to WFPC2 as in Fig. 1131 to those from fitting models to the galaxies. As we see in Fig. 1131 the luminosity-weighted moment 
estimates are generally rounder than the ones measured by GALFIT model fitting, despite the higher resolution and better sampling of 
ACS as compared with WFPC2. This figure shows that modeling the PSF when measuring the eUipticity is particularly important for 
deriving a robust estimate of t he u nderlying value, and that the resolution of the instrument is a lesser factor (at least in this case between 
WFPC2 and ACS). As in Fig. 1131 the galaxies classified as SOs are more likel y to have small eUipticities when they are larger in size. The 
median offset for SO galaxies is <5emed = —0.10, smaller than we found in Fig. 1131 but still much larger than the statistical and systematic 
errors we find for our PSF-corrected GALFIT model approach. 

In FigurefMl we plot the eUipticities we measure using SExtractor on the ACS images of MS 1054-03 and CL 0024-1-16 
in comparison to the eUipticities we find by fitting models using GALFIT. From these plots, we conclude that the 
more compact PSF for ACS is not the dominate reason for the different eUipticities between our work and the work 
of previous authors. Removing the effect of the PSF is the most important component. 

Finally, we plot the comparison between the eUipticities as derived by GALFIT and those from SExtractor in Figure 
[T5l for our simulations of high- redshift galaxies. These are images of real galaxies that have been redshifted and 
to which we have added noise. As with the previous plots showing trends in the data, the simulated galaxies are 
measured to be rounder with SExtractor than by fitting models using GALFIT. The average offset is Semed ~ —0.09 
between the SExtractor eUipticity and the GALFIT eUipticity, in good agreement with the Scmed — —0.10 we measure 
in Figure [TH As we show in Figure [H our GALFIT eUipticities reproduces the underlying eUipticity distribution. 
Because our GALFIT method reproduces the input eUipticities so well, we conclude that the SExtractor method of 
measuring eUipticities underestimates the eUipticities of galaxies at higher redshifts. We would advise against the use 
of SExtractor eUipticities for quantitative studies of high redshift galaxies using ACS and WFPC2 where the PSF is 
a modest, but significant fraction of the size of the object. 

B. THE REDSHIFT CATALOG FOR CL J1226+33 

Of the clusters discussed h ere, only CL J1226-I- 33 has no previously published redshift catalog beyond the six galaxies 
listed in the discovery paper (|Ebeling et al.|[200ll ). We have collected redshifts for 21 new members in the central region 
of the cluster, along with 10 more galaxies outside of the region covered by the ACS imaging. 30 of these 31 galaxies 
were selected to have ACS Ign magnitudes < 24 mag AB and colors in the range 1.5 < Veoe — < 2, the remaining 
one was a "filler" object with a bluer color. The data were taken using the DEIMOS spectrograph on Keck II, in a 1 
hr exposure with the 600 line mm~^ grating in variable conditions. We list the catalog of redshifts in Table |4l 

We measured redshifts by centroiding emission features and cross-correlating with a variety of templates. We 
computed the biweight center and scale of all galaxies with redshifts within 0.87 < z < 0.91 . We found z = . 890ib 0.001 
with a dispersion of 1143 ± 162 km s~^, in good agreement with the 1270 km s^^ quoted in |j0rgensen et al] (|2006l ). We 
used a biweight to estimate both the redshift and the velocity dispersion of the cluster with errors from the jackknife 
of the galaxy redshifts. The redshift range we used for membership excludes two galaxies at 2 ~ 0.91. Including those 
two objects raises the dispersion to 1322 ± 221 km s^^. It is often the case that, with small numbers of galaxies 
that the dispersion is overestimated - see for example the dispersions quoted for CL 1604-1-4304 and CL 1604-f4321 
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Fig. 15. — SExtractor ellipticities as a function of the GALFIT ellipticities for galaxies from our simulations. We plot the SExtractor 
ellipticities on the y-axis as a function of the GALFIT ellipticities fo r th e same simulated galaxies. For comparison, we plot a line with 
a slope of one. Note that this plot shows a similar behavior as Figs. 1131 and 1141 namely that the ellipticities as estimated by SExtractor 
are rounder than those from GALFIT. The GALFIT measurements of the simulated z > 0.3 galaxies recover the actual ellipticities with 
an average offset of ht ~ —0.01 (see Fig. [TJ, much smaller than the offset of ~ —0.09 we find between the GALFIT and SExtractor 
ellipticities. The higher density of points in the z = 0.85 redshift bin simply reflects the larger number of clusters included in that bin. 

m iLubin et al.l (|2000f ) versus those from the much larger sample in iGal et all (|2005D . This arises because nearby large 
scale structure or groups of galaxies can project into line of sight of the cluster. We opted to remove the two galaxies 
at 2; ~ 0.91 to provide conservative estimate of the dispersion, which we also list in Tabled 
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TABLE 4 
Catalog for CL J1226+33 



Obj 


R.A. 


Dec. 


I 


V-I 


z 




(J2000) 


(J2000) 


(mag AB) 


mag AB 





GAL12263974- 


h3329315 


12 


:26: 


39 


.74 


33 


:29 


:31 


.51 


21 


.82 


± 


0, 


,02 


1 


.75 


± 


0, 


,03 


1 


.035 


GAL12264240- 


h3330094 


12 


:26: 


42 


.40 


33 


:30 


:09 


.40 


22 


.11 


± 


0, 


,01 


1 


.67 


± 


0, 


,01 





.929 


GAL12263338- 


h3330277 


12 


:26: 


33 


.38 


33 


:30 


:27, 


.77 


21 


.68 


± 


0, 


.01 


1 


.82 


± 


0, 


,02 





.911 


GAL12264935- 


h3330312 


12 


:26: 


49 


.35 


33: 


:30: 


:31 


.25 


21 


.61 


± 


0, 


,01 


1 


.82 


± 


0, 


,01 


0, 


.891 


GAL12263934- 


h3330353 


12 


:26: 


:39 


.34 


33: 


:30: 


:35, 


.35 


22 


.40 


± 


0, 


,01 


1 


.96 


± 


0, 


,03 


0, 


.943 


GAL12264600- 


h3330347 


12 


:26: 


:46 


.00 


33: 


:30: 


:34, 


.74 


21 


.43 


d= 


0, 


,01 


1 


.84 


± 


0, 


,01 


0, 


.892 


GAL12263157- 


h3330473 


12 


:26: 


31 


.57 


33 


:30 


:47, 


.38 


22 


.99 


± 


0, 


,02 


0, 


.72 


± 


0, 


,02 


0, 


.963 


GAL12262950- 


h3330580 


12 


:26: 


29 


.50 


33 


:30 


:58 


.01 


21 


.59 


± 


0, 


,01 


1 


.81 


± 


0, 


,02 


1 


.211 


GAL12264172- 


h3331041 


12 


:26: 


41 


.72 


33: 


:31 


:04, 


.15 


22, 


.26 


± 


0, 


,01 


0, 


.59 


± 


0, 


.01 


0, 


.690 


GAL12264671- 


h3331037 


12 


:26: 


:46, 


.71 


33: 


:31 


:03, 


.73 


22, 


.07 


± 


0, 


.01 


1 


.60 


± 


0, 


.01 


0, 


.767 


GAL12263998+3331047 


12 


;26; 


;39, 


.98 


33: 


;31 


;04, 


.72 


21, 


.70 


± 


0, 


.01 


1 


.85 




0, 


.02 


0, 


.929 


GAL12263648+3331124 


12: 


:26; 


;36 


.48 


33: 


:31 


:12, 


.47 


22, 


.47 


lb 


0, 


.02 


0, 


.99 


± 


0, 


.02 


0, 


.846 


GAL12263107+3331349 


12: 


;26: 


:31, 


.07 


33: 


;31: 


;34, 


.97 


22, 


.73 


± 


0. 


.01 


0, 


.60 


± 


0. 


.01 


0, 


.965 


GAL12263402+3331356 


12: 


;26: 


:34, 


.02 


33: 


;31: 


;35, 


.62 


21, 


.84 


± 


0. 


.01 


1, 


.84 


± 


0. 


.02 


0, 


.966 


GAL12265120+3331385 


12: 


:26: 


:51, 


.20 


33: 


;31: 


;38. 


.54 


22, 


.68 


± 


0. 


.01 


1, 


.90 


± 


0. 


.02 


0, 


.881 


GAL12265293+3331465 


12: 


;26: 


:52, 


.93 


33: 


;31: 


:46. 


.54 


22. 


.78 


± 


0. 


.01 


1, 


.76 


± 


0. 


.02 


0, 


.964 


GAL12263832+3331482 


12: 


:26: 


:38, 


.32 


33: 


;31: 


:48. 


.21 


23. 


.27 


± 


0. 


.06 





.63 


± 


0. 


.09 


0, 


.542 


GAL12262701+3331580 


12: 


:26; 


:27, 


.01 


33: 


:31: 


:58, 


.09 


22, 


.27 


± 


0, 


.08 


1, 


.82 


± 


0, 


.08 


0, 


.645 


GAL12272134- 


h3332083 


12 


:27: 


:21 


.34 


33: 


:32: 


:08, 


.32 


21 


.86 


± 


0, 


,01 


1 


.81 


± 


0, 


,01 


0, 


.893 


GAL12262817- 


h3332232 


12 


:26: 


:28 


.17 


33 


:32 


:23 


.23 


22 


.93 


± 


0, 


,02 


1 


.13 


± 


0, 


,03 





.673 


GAL12271941- 


h3332269 


12 


:27: 


:19 


.41 


33 


:32: 


:26, 


.91 


23 


.08 


± 


0, 


,02 


0, 


.70 


± 


0, 


,02 





.885 


GAL12265525- 


h3332324 


12 


:26: 


:55, 


.25 


33: 


:32: 


:32, 


.43 


22 


.28 


± 


0, 


,01 


1 


.90 


± 


0, 


,01 


0, 


.895 


GAL12271460- 


h3332377 


12: 


:27: 


14, 


.60 


33: 


:32: 


:37, 


.70 


22, 


.10 


± 


0, 


,01 


1 


.89 


± 


0, 


,01 


0, 


.881 


GAL12265629- 


h3332414 


12 


:26: 


56, 


.29 


33: 


:32: 


:41 


.48 


22 


.53 


± 


0, 


,01 


1 


.78 


± 


0, 


,01 


0, 


.893 


GAL12265995- 


h3332405 


12 


:26: 


59, 


.95 


33: 


:32: 


:40, 


.54 


22 


.48 


± 


0, 


,01 


1 


.76 


± 


0, 


,02 


0, 


.892 


GAL12265923- 


h3332405 


12 


:26: 


:59 


.23 


33 


:32 


:40 


,59 


22 


.21 


± 


0, 


,01 


1 


.82 


± 


0, 


,01 


0, 


.897 


GAL12265689- 


h3332437 


12 


:26: 


56 


.89 


33 


:32 


:43 


,75 


22 


.86 


± 


0, 


,01 


1 


.84 


± 


0, 


,02 





.896 


GAL12270510- 


h3332475 


12 


:27: 


05 


.10 


33 


:32 


:47, 


.53 


23 


.01 


± 


0, 


,01 





.77 


± 


0, 


,01 





.798 


GAL12265060- 


h3332461 


12 


:26: 


:50, 


.60 


33: 


:32: 


:46, 


.18 


21, 


.50 


± 


0, 


.01 


1 


.83 


± 


0, 


.01 


0, 


.875 


GAL12265825+3332485 


12: 


;26: 


:58, 


.25 


33: 


;32: 


:48, 


.57 


19, 


.13 


± 


0. 


.00 


1, 


.96 


± 


0. 


.00 


0, 


.891 


GAL12271547+3332539 


12: 


;27: 


:15, 


.47 


33: 


;32: 


:53, 


.91 


22, 


.87 


± 


0. 


.02 


0, 


.69 


± 


0. 


.02 


1, 


.034 


GAL12270083+3333019 


12: 


;27: 


:00, 


.83 


33: 


:33: 


;01. 


.97 


22, 


.50 


± 


0. 


.01 


1, 


.35 


± 


0. 


.02 


0, 


.914 


GAL12265214+3333071 


12: 


;26: 


:52, 


.14 


33: 


;33: 


:07. 


.12 


23, 


.02 


± 


0. 


.02 


1, 


.76 


± 


0. 


.02 


0, 


.881 


GAL12265714+3333046 


12: 


;26: 


:57, 


.14 


33: 


;33: 


;04, 


.61 


21, 


.65 


± 


0. 


.01 


1, 
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